Interaction-Induced Strong Localization in Quantum Dots 
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We argue that Coulomb blockade phenomena are a useful probe of the cross-over to strong correla- 
tion in quantum dots. Through calculations at low density using variational and diffusion quantum 
Monte Carlo (up to r s ~ 55), we find that the addition energy shows a clear progression from fea- 
tures associated with shell structure to those caused by commensurability of a Wigner crystal. This 
cross-over (which occurs near r s ~ 20 for spin-polarized electrons) is, then, a signature of interaction- 
driven localization. As the addition energy is directly measurable in Coulomb blockade conductance 
experiments, this provides a direct probe of localization in the low density electron gas. 



Localization of electrons induced by electron-electron 
interactions is a key issue in strongly-interacting systems 
[l[ which dates from the dawn of solid state physics when 
Wigner introduced the notion of an electron crystal Q. 
The recent focus on localization in inhomogeneous sys- 
tems is stimulated in part by experiments which appar- 
ently see a metal-insulator transition in two dimensions 
This occurs in low-density electron gas samples in 
which disorder is presumably important. Recall that in a 
degenerate electron gas, low density implies strong inter- 
actions. Here we show that an approach from nanoscale 
physics, in which the potential confining the electrons 
plays a key role analogous to that of disorder, provides 
new information on this critical problem. 

The Coulomb blockade effect has been a valuable tool 
for probing a variety of interaction effects in quantum 
dots [1, [1, d, 0- By "quantum dot" we mean a con- 
fined region of electron gas containing between N = 1 
and ~ 1000 electrons; experimentally, they have proven 
remarkably tunable through the use of gate voltages on 
nearby electrodes. We treat N < 20 here. The large elec- 
trostatic charging energy usually forces the dot to have 
a fixed number of electrons, preventing the flow of cur- 
rent. The blockade is lifted by using a gate to tune the 
energy for N electrons to be the same as that for N + 1 , 
inducing a finite conductance. The conductance through 
the dot as a function of gate voltage is therefore a series 
of sharp peaks. The height and position of these peaks 
encode information about the dot's ground state; for in- 
stance, the spacing between the peaks is proportional to 
the second difference of the ground-state energy with re- 
spect to electron number N, a quantity known as the 
addition energy. Quantum many-body physics probed in 
quantum dots includes, for instance, the atomic-like ef- 
fect of exchange and correlation in altering the filling of 
single-particle shells . several kinds of Kondo effects 
, aspects of the fractional quantum Hall effect [H , and 
the entanglement of spin and orbital degrees of freedom 
for quantum information purposes 

We show that the Coulomb blockade can be a valuable 
probe of localization in quantum dots. In particular, the 



characteristic pattern of the addition energy changes as 
the dot crosses over from the extended states of the high 
density Fermi liquid regime to the localized states char- 
acterizing the low density Wigner crystal. Interaction 
strength is often characterized by r s = 1 /a* B \J im (in two 
dimensions) where a* B is the effective Bohr radius and n is 
the electron density. In our quantum dots, the cross-over 
occurs for r s substantially smaller than the value [slli^ lloj 
in bulk. This decrease in interaction strength needed for 
localization is connected to the density inhomogeneity 
necessarily present in this confined system; density in- 
homogeneity produced in other ways, such as by disor- 
der, may similarly enhance localization. Our main point 
is that such a cross-over to localization can be directly 
measured in Coulomb blockade experiments. 

The cross-over from Fermi liquid to "Wigner molecule" 
in quantum dots was studied previously using various 
many-body methods. Exact diagonalization 0,[fl|, while 
the most robust and direct approach, is limited to small 
electron number and small r s due to conver gen ce prob- 
lems. Path- integral quantum Monte Carlo 12j, |l3, [l4| 
is well suited for finite temperature properties, but pre- 
serves only S z symmetry and has large statistical fluc- 
tuations at low temperatures. Variational and diffusion 
Monte Carlo [l5| preserve S z , S 2 , and L z , symmetry; 
though limited by the "fixed node" error which depends 
on the quality of the trial wave function [l5j], they were 
successfully used to study quantum dots for r s up to ~ 4 



[liL Hyl IS. [l9l| and then up to r s ~ 15 in our recent 
work [20j, l2l| . Here, we apply recently developed energy 



minimization methods [22j, [23j to floating gaussian-based 



trial wave functions, enabling us to decrease the fixed- 
node error and investigate the strongly correlated regime 
up to as high as r s ^55. 

Our model quantum dot consists of N interacting elec- 
trons in a two-dimensional circular quadratic potential. 
The Hamiltonian, expressed in effective atomic units 
(electronic charge e, dielectric constant e, effective mass 
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m*, and h are set to 1), is given by 
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where u> is the spring constant of the quadratic potential 
which provides control of the strength of the Coulomb 
interaction with respect to the kinetic energy. In analogy 
with 2D bulk systems, we characterize the interaction 
strength using r s — (ttH)^ 1 ^ 2 , where n = J n 2 (r)dr/N is 
the mean density of electrons. 

Variational (VMC) and diffusion (DMC) Monte Carlo 
techniques [l5| were used to calculate the properties of 
our model quantum dots. One starts with a set of single- 
particle orbitals — simple gaussian functions or from a 
self-consistent calculation (Hartree or Kohn-Sham). We 
then perform a VMC calculation using a trial wave func- 
tion, which is a linear combination of products of 
up- and down- spin Slater determinants of these orbitals 
multiplied by a Jastrow factor. (The detailed form of our 
Jastrow factor is in Ref. [H.) In a second stage, we use 
fixed-node DMC 15|, l24| to project the optimized many- 
body wave function onto a better approximation of the 
true ground state, an approximation that has the same 
nodes as $>t- The fixed- node DMC energy is an upper 
bound to the true energy and depends only on the nodes 
of the trial wave function obtained from VMC. 

In order to capture the Wigner molecule aspect, we 
build our trial wave functions from floating gaussian or- 
bitals. The use of localized orbitals in confined structures 
causes three problems: (i) The positions of the gaussians 
are initially unknown and must be optimized together 
with the widths, (ii) Slater determinants built from lo- 
calized gaussians do not conserve angular momentum un- 
less integrated over all angles j2^|. Strikingly, we find 
that taking a linear superposition of as few as 3 rotated 
determinants (for each ring) recovers more than 99% of 
the missing DMC energy, (iii) Low spin states require 
the superposition of many Slater determinants in order 
to have an exact spin eigenstate, making spin polarized 
states easier to study. All the parameters in our VMC 
calculations — Jastrow parameters, widths and positions 
of the gaussians, and determinantal coefficients — are op- 
timized using efficient energy minimization techniques 
[22I [23| . The use of these techniques is a key feature 
which ultimately allows us to cover a much wider range 
of density with this method than was previously possible. 

We have also used liquid-like orbitals obtained by solv- 
ing Kohn-Sham or Hartree equations. A comparison of 
the results obtained from different types of orbitals is 
presented below. Quantities that do not commute with 
the Hamiltonian were calculated u sing an extrapolated 
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estimator Fqmc = 2Fdmc - Fvmc 

In Fig. Q] we plot the addition energy A 2 E(N) = 
E G (N + 1) + E G (N - 1) - 2E G (N) of the fully spin po- 
larized ground states for different lj values, including the 
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FIG. 1: Addition energy as a function of number of electrons 
in the quantum dot for different ui values (correspondin g Vs 
values are approximately 0, 2, 5, 15, 25, 55). The classical |2q ] 
and non-interacting addition energies are also shown. Note 
the progression from peaks determined by shell structure to 
those consistent with classical electrostatics. 



results for the classical [26j and non- interacting systems. 
As the addition energy is directly measurable in Coulomb 
blockade conductance experiments [H, 0|, the structure 
of the curves in Fig. [1] can be observed by applying a 
strong in-plane magnetic field. In the non-interacting 
case (uj — > 00), the peaks at N = 3, 6, 10, and 15 are due 
to the shell structure of the two-dimensional harmonic 
confinement potential. As u> is decreased, the electronic 
density decreases and electron-electron interactions be- 
come more important. Note that this causes the strength 
of the fluctuations in A 2 E to decrease. Around uj ~ 0.015 
(r s ~ 15), the shell structure peaks are washed out by the 
Coulomb interaction energy. For larger r s values, a new 
structure sets in, showing peaks at N = 3, 7, 11, 13, 15, 
and 18. These new magic numbers are caused by com- 
mensurability of a Wigner crystal. In fact, at r s ~ 55, 
there is an impressive quantitative agreement between 
the quantum mechanical and classical addition energies. 
Our results suggest that the cross-over from a liquid-like 
state towards the Wigner crystal regime in quantum dots, 
which occurs around r s ~ 20 for fully polarized electrons, 
can be directly observed experimentally. 

In the 2D electron gas, QMC methods were previously 
used to study the liquid-solid transition by comparing 
the fixed-node DMC energies as a function of r s obtained 
using liquid-like orbitals and localized gaussian orbitals 
[S @> E3] • The energies using the two types of orbitals 
cross near r s ~ 35, which was therefore identified as the 
localization point. In order to investigate if such a cross- 
over occurs in quantum dots, we compared the DMC 
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FIG. 2: Difference between the DMC energies obtained us- 
ing liquid-like (Hartree) orbitals and solid-like (localized gaus- 
sian) orbitals as a function of r a . Both types of orbitals appear 
to work well throughout the whole range of r s . 



energies obtained using delocalized Hartree orbitals and 
localized gaussian orbitals for several (N, S) states. Fig. [5] 
shows the energy differences for 4 different cases. To our 
surprise, we did not find any systematic r s dependence: 
solid and liquid type orbitals give equally good DMC 
energies at r s values ranging from ~2 to 55. 

We now consider an alternative way of quantifying the 
degree of localization as a function of r s : We define the 
angular power spectrum of the quantum dot as 
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where r = (r, 9) is the position vector and Tg is the 
Fourier transform in the angular direction. For a classical 
ring of N r equally spaced electrons, the above expression 
is proportional to Ylm=i i^r — m ) cos(2itmkg /N r ), which 
has maxima at kg = 0, N r , etc. 

Fig. [3[a) shows the power spectrum of a fully polar- 
ized system of N — 18 electrons at r s — 52. From 
this plot, three important quantitative points can be ex- 
tracted: (i) At kg = 0, Eq. ([2]) reduces to the probability 
of finding two electrons at a given radius; thus, the peaks 
at kg — give the radial location of the electrons. There 
is an inner ring at r ~ 110 well separated from the outer 
ring at r ~ 230 (plus a single electron at the center), (ii) 
There are two more peaks at kg — 6 and kg = 11, indicat- 
ing the localization of 6 and 11 electrons in the inner and 
outer rings, respectively, (iii) The height of the secondary 




FIG. 3: Angular power spectrum for N — 18 as a function of 
ke and quantum dot radius r for (a) r s — 52 and (b) r 3 = 4.8. 
Insets: Corresponding pair-densities (fixed electron marked 
by "x"). From the peaks in the power spectrum, one deduces 
the number of electrons in each ring and the strength of the 
angular modulation. 



peaks compared to the kg = peaks is a measure of the 
angular localization of the electrons in the rings. For 
comparison, the pair-density is shown in the inset — the 
density of electrons at r given that one electron is fixed 
at the position marked with the cross. We note that the 
(1,6,11) configuration (that is, 1 electron at the center, 
6 in the inner ring, and 11 in the outer ring) indicated 
by both the power spectrum and pair-density plots also 
corresponds to the classical configuration [26j. Clearly, 
at r s ~ 52, the electrons are well localized at their classi- 
cal position, confirming the agreement between quantum 
mechanical and classical addition energies in Fig. [1] 

On the other hand, Fig.[3jb) shows the results for r s ss 
5. The outer and inner rings are not separated as well, 
and the peaks at kg = 6 and kg = 11 are much weaker, 
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FIG. 4: The height of the peaks in the angular power spec- 
trum at r = r max for 6 electrons in three cases: fully polar- 
ized (S = 3), unpolarized, and S = 2 partially polarized. For 

5 = 0, the number of electrons in the outer ring changes from 

6 to 5 as the strength of the interactions increases. 



although still visible despite the relatively low r s . The 
pair-density shows weak intrashell angular modulation 
consistent with the kg = 11 peak and almost no intershell 
angular correlation. 

For small particle number, N < 8, we studied several 
spin states of the quantum dot, not just the fully po- 
larized case. N — 6 is particularly interesting: two spa- 
tial configurations of the electrons have similar energies — 
there are either 6 electrons in one ring or 5 electrons in 
a ring around 1 electron. Fig. [4] shows the development 
of angular modulation in this case. For full polariza- 
tion, the amplitude of the (1,5) configuration smoothly 
increases. We find this type of smooth increase for all 
values of N when fully polarized. In contrast, for 5 = 0, 
(0,6) is dominant at small r s while (1,5) takes over for 
strong interaction. For r s > 20, the amplitude of the 
modulation is identical in the two cases, as well as for 
the intermediate S — 2 case. This intriguing insensitiv- 
ity of the spatial behavior to the spin state is another 
indication of localization. 

In conclusion, we have proposed a new way of prob- 
ing Wigner localization in quantum dots and presented 
three kinds of evidence for a cross-over from extended to 
localized electronic states at around r s ~ 20. First, the 
angular power spectrum develops sharp peaks indicating 
that the positions of the electrons become highly corre- 
lated (Fig. [3]). Second, for N = 6 in which two spatial 
configurations of the electrons compete, the angular mod- 
ulation becomes completely independent of the spin state 
(Fig. HJ. Finally and most importantly, structure in the 
addition energy for spin-polarized electrons as a function 
of electron number shows a clear progression from peaks 
consistent with shell structure to those consistent with 



the electrostatics of point particles (Fig. [TJ. This latter 
effect should be observable in Coulomb blockade conduc- 
tance measurements, allowing for the first time a direct 
probe of electronic localization in a two-dimensional sys- 
tem. 

We thank D. Ullmo and V. Elser for valuable conver- 
sations. This work was supported in part by the NSF 
(DMR-0506953) and by the DOE-CMSN (DE-FG02- 
07ER46365) grant. AG was supported in part by funds 
from the David Saxon Chair at UCLA. 



[1] 

[2 
[3 

[4 

[5 

[6 

ir 

[8 
[9 

[10 
[11 

[12 

[13 

[14, 

[is; 

[16 

[i7; 
[is: 

[19 
[20 
[21 
[22 
[23 



S. A. Kivelson and E. Fradkin, in Handbook of High- 
Temperature Superconductivity, edited by J. R. Schri- 
effer and J. S. Brooks (Springer- Verlag, Berlin, 2007), 
arXiv:cond-mat/0507459. 

E. Wigner, Phys. Rev. 46, 1002 (1934). 

S. V. Kravchenko and M. P. Sarachik, Rep. Prog. Phys. 
67, 1 (2004). 

D. Heiss, ed., Quantum Dots: A Doorway to Nanoscale 
Phystcs (Springer- Verlag, Berlin, 2005). 
S. M. Reimann and M. Manninen, Rev. Mod. Phys. 74, 
1283 (2002), and references therein. 

L. P. Kouwenhoven, D. G. Austing, and S. Tarucha, Rep. 
Prog. Phys. 64, 701 (2001). 

M. Gravis, I. G. Rau, R. M. Potok, and D. Goldhaber- 
Gordon, in Handbook of Magnetism and Advanced Mag- 
netic Materials, edited by H. Kronmiiller and S. Parkin 
(Wiley, New York, 2007). 

B. Tanatar and D. M. Ceperley, Phys. Rev. B 39, 5005 
(1989). 

C. Attaccalite, S. Moroni, P. Gori-Giorgi, and G. B. 
Bachelet, Phys. Rev. Lett. 88, 256601 (2002), and 91, 
109902(E) (2003). 

X. Waintal, Phys. Rev. B 73, 075417 (2006). 

M. Rontani, C. Cavazzoni, D. Bellucci, and G. Goldoni, 

J. Chem. Phys. 124, 124102 (2006). 

R. Egger, W. Hausler, C. H. Mak, and H. Grabert, Phys. 
Rev. Lett. 82, 3320 (1999), and 83, 462(E) (1999). 
A. V. Filinov, M. Bonitz, and Y. E. Lozovik, Phys. Rev. 
Lett. 86, 3851 (2001). 

S. Weiss and R. Egger, Phys. Rev. B 72, 245301 (2005). 
W. M. C. Foulkes, L. Mitas, R. J. Needs, and G. Ra- 
jagopal, Rev. Mod. Phys. 73, 33 (2001). 

F. Pederiva, C. J. Umrigar, and E. Lipparini, Phys. Rev. 
B 62, 8120 (2000), and 68, 089901(E) (2003). 

A. Harju, S. Siljamaki, and R. M. Nieminen, Phys. Rev. 
B 65, 075309 (2002). 

A. D. Guglii, G. S. Jeon, C. J. Umrigar, and J. K. Jain, 
Phys. Rev. B 72, 205327 (2005). 

A. Ghosal, C. J. Umrigar, H. Jiang, D. Ullmo, and H. U. 

Baranger, Phys. Rev. B 71, 241306(R) (2005). 

A. Ghosal, A. D. Guglii, C. J. Umrigar, D. Ullmo, and 

H. U. Baranger, Nature Phys. 2, 336 (2006). 

A. Ghosal, A. D. Giiglii, C. J. Umrigar, D. Ullmo, and 

H. U. Baranger, Phys. Rev. B (2007), in press. 

C. J. Umrigar and C. Filippi, Phys. Rev. Lett. 94, 150201 

(2005). 

C. J. Umrigar, J. Toulouse, C. Filippi, S. Sorella, and 
R. G. Hennig, Phys. Rev. Lett. 98, 110201 (2007). 



■5 



[24] C. J. Umrigar, M. P. Nightingale, and K. J. Runge, J. [26] V. M. Bedanov and F. M. Peeters, Phys. Rev. B 49, 2667 

Chem. Phys. 99, 2865 (1993). (1994). 
[25] Y. S. Li, C. Yannouleas, and U. Landman, Phys. Rev. B 

73, 075301 (2006). 



